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The irreversible growth of a binary mixture under far-from-equilibrium conditions is studied in 
three-dimensional confined geometries of size Lx x Ly x Lz, where Lz ^ Lx = Ly is the growing 
direction. A competing situation where two opposite surfaces prefer different species of the mixture 
is analyzed. Due to this antisymmetric condition an interface between the different species develops 
along the growing direction. Such interface undergoes a localization-delocalization transition that 
is the precursor of a wetting transition in the thermodynamic limit. Furthermore, the growing 
interface also undergoes a concave-convex transition in the growth mode. So, the system exhibits 
a multicritical wetting point where the growing interface is almost flat and the interface between 
species is essentially localized at the center of the film. 



The study of interfacial phenomena of materials con- 
fined in thin film geometries under equilibrium condi- 
tions has drawn enormous attention over the last decades 
[1-3]. Interfacial phase transitions and critical phenom- 
ena in confined samples exhibit a quite distinct physical 
behavior compared to that occurring in the bulk due to 
the finite separation between the walls and the specific 
wall-particle interaction [1]. The understanding of this 
new type of phenomena is not only of general and fun- 
damental interest but may also play an important role 
in the development of new technologies. Accordingly, 
the properties of films of magnetic materials [4-7] , poly- 
mer blends [8,9], binary mixtures [10], fluids [11], etc, 
have been investigated thoroughly. Liquid-gas condensa- 
tion under confinement between identical walls displays 
capillary condensation when the walls promote the fluid 
phase [12]. Mapping the fluid into a lattice gas model 
and using the magnetic terminology, the wall- particle 
interaction is accounted for by a surface magnetic field 
[4-7]. For this reason, Ising-like models are very useful 
to understand the underlying physics of more complex 
systems [13]. Another interesting scenario takes place 
when an Ising film is confined between two competing 
walls a distance L apart from each other, so that the 
surface magnetic fields (H) are of the same magnitude 
but opposite direction. These competing fields cause the 
emergence of an interface that undergoes a localization- 
delocalization transition at an L— dependent temperature 
Tu,(L, H) that is the precursor of the true wetting tran- 
sition temperature Tw{H) of the infinite system. 

It is surprising that, in contrast to the effort devoted 
to the understanding of wetting phenomena under equi- 
librium conditions, the non-equilibrium counterpart has 
received very little attention. For instance, Hinrichsen et 
al. [14] have introduced a nonequilibrium growth model 
of a one dimensional interface that undergoes a transition 
from a bounded state to a non-bounded one. 

Within this context, the aim of this work is to perform 
an extensive numerical study of the irreversible growth 
of a magnetic material confined between parallel walls 
where competing surface magnetic fields act. It is shown 



that the interplay between confinement and growth mode 
leads to a physically rich phase diagram (on the plane 
H — T) that exhibits a localization-delocalization tran- 
sition in the interface that runs along the walls and a 
change of the curvature of the growing interface running 
perpendicularly to the walls. Extrapolation of this sce- 
nario to the thermodynamic limit leads to a multicritical 
wetting point. 

The growth of a ferromagnetic material, with spins 
having two possible orientations, is studied using the 
so called Magnetic Eden Model (MEM) [15]. Monte 
Carlo simulations are performed on the square lattice 
in (2 + 1)— dimensions, using a rectangular geometry 
Lx X Ly X Lz with L^ ^ L^ = Ly = L. The location of 
each site on the lattice is specified through its rectangu- 
lar coordinates {i,j,k), (1 < i,j < L, 1 < k < L^)- The 
starting seed for the growing cluster is a plane oi L x L 
spins placed at A: = 1 and cluster growth takes place along 
the positive longitudinal direction (i.e.. A: > 2). The 
boundary conditions are periodic along one of the trans- 
verse directions (say in the i— direction) but open along 
the remaining transverse direction. In the latter, compet- 
ing surface magnetic fields Hi > {LLl < 0) acting on 
the sites placed at j = 1 {j — L), with H = Hi = \Hl\, 
are considered. Then, clusters are grown by selectively 
adding spins {Stjk = ±1) to perimeter sites, which are 
defined as the nearest-neighbor (NN) empty sites of the 
already occupied ones. 

Considering a ferromagnetic interaction of strength 
J > between NN spins, the energy E oi a, given config- 
uration of spins is given by 
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where the summation {ijk, i j k ) is taken over occupied 
NN sites, while {ik, Si), {ik^Yji) denote summations car- 
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ried over occupied sites on the surfaces j = 1 and j = L, 
respectively. Throughout this work we set the Boltz- 
mann constant equal to unity and we take the temper- 
ature, energy, and magnetic fields measured in units of 
J. The probability for a perimeter site to be occupied 
by a spin is taken to be proportional to the Boltzmann 
factor exp(— -^P), where AE is the change of energy in- 
volved in the addition of the given spin. At each step, the 
probabilities of adding up and down spins to a given site 
have to be evaluated for all perimeter sites. After proper 
normalization of the probabilities, the growing site and 
the orientation of the spin are determined with Monte 
Carlo techniques. Using this procedure, clusters having 
up to 10^ spins have typically been grown. Although 
both the interaction energy and the Boltzmann proba- 
bility distribution considered for the MEM are similar 
to those used for the Ising model with surface magnetic 
fields [9,12], it must be stressed that these two models 
operate under extremely different conditions, namely the 
MEM describes the irreversible growth of a magnetic 
material and the Ising model deals with a magnet under 
equilibrium. Previous studies have demonstrated that 
the MEM in (1 -|- 1)— dimensions is not critical but it ex- 
hibits a second-order transition at Tc = 0.69 ± 0.01 in 
(2 + 1)— dimensions [16]. 

The phase diagram of the MEM in a confined geome- 
try with competing surface fields is very rich and exhibits 
eight regions. For the sake of clarity we will first discuss 
the main characteristics of each of these regions by means 
of snapshot configurations (see figure 1). Subsequently, 
we will quantitatively locate the boundary between these 
regions and plot the corresponding phase diagram (see 
figure 2). Finally, we will draw the phase diagram in the 
thermodynamic limit (see inset of figure 2) by extrapo- 
lating finite-size results (see figure 3). 

The vertical straight line at Tc{L = 12) = 0.84 (see 
figure 2) represents the L— dependent "critical" temper- 
ature of the finite system [17]. So, the left (right) hand 
side part of the phase diagram corresponds to the or- 
dered (disordered) growth regime that involves Regions 
/, //, ///, IV, and A (Regions V, VI, and B). 

For low temperatures and small fields (Region / of fig- 
ure 2, see also figure 1(a)), the system grows in an ordered 
state and a domain having mostly spins with a single 
orientation prevails. Due to thermal fiuctuations, small 
clusters with the opposite orientation may appear, prefer- 
ably on the surface where the field competing with the 
dominant orientation is applied. These "drops" might 
grow and when the typical size of the fluctuation is of 
the order of L, a magnetization reversal may occur, thus 
changing the sign of the dominant domain. The forma- 
tion of sequences of well ordered domains is a finite size 
effect that is also characteristic of the ordered phase of 
some confined spin systems such as the Ising magnet [7]. 
Due to the open boundary conditions, empty perimeter 
sites at the confinement walls experience a missing neigh- 
bor effect. Since H is too weak to compensate this ef- 
fect, the system grows preferentially along the center of 



the sample as compared to the walls, and the resulting 
growth interface exhibits a convex shape. So, Region / 
corresponds to the Ising-like nonwet state and the convex 
growth regime. 
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FIG. 1. Snapshot pictures showing a longitudinal slice 
given by a fixed value of the transverse coordinate i. Red 
(black) circles correspond to spins up (down). The surface 
field on the upper (lower) confinement wall is positive (neg- 
ative). The snapshots correspond to a lattice size L = 32 
and several different values of temperature and surface fields: 
{a.)H = 0.05, T = 0.6; {h)H = 0.5, T = 0.55; {c)H = 1.4, 
T = 0.6; {d)H = 0.1, T = 1.0; {e)H = 1.6, T = 1.4; and (f) 
H = 0.20, T = 0.82. 

Increasing the field but keeping the temperature low 
enough, one may enter Region // (figure 2, see also fig- 
ure 1(b)). Here the system is still in its ordered phase and 
neighboring spins grow preferably parallel-oriented. The 
surface fields in this region are stronger and thus capa- 
ble of compensating the missing NN sites at the surfaces. 
But, since the fields on both surfaces have opposite ori- 
entations, one has that, on the one hand, the field that 
has the same orientation as that of the dominant spin 
cluster will favor the growth of surface spins, while on 
the other hand, the sites on the surface with opposite 
field will have a lower growing probability. Thus, on this 
disfavored side the growing interface becomes pinned and 
the curvature of the growing interface is not defined. 

Keeping H fixed within Region // but increasing the 
temperature, thermal noise will enable the formation of 
drops on the disfavored side that eventually may nucle- 
ate into larger clusters as the temperature is increased 
even further. This process may lead to the emergence 
of an up-down interface, separating oppositely oriented 
domains, running in the direction parallel to the walls. 
Since sites along the up-down interface are surrounded 
by oppositely oriented NN spins, they have a low growing 
probability. So, in this case the system grows preferably 
along the confinement walls and the growing interface 
is concave (figure 1(c)). Then, as the temperature is in- 
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creased, the system erosses to Region A (see figure 2) and 
we observe the onset of two competitive growth regimes: 
(i) one exhibiting a non-defined growing curvature that 
appears when a dominant spin orientation is present, as 
in the case shown in figure 1(b); (ii) another that appears 
when an up-down interface is estabhshed and the system 
has a concave growth interface, as is shown in figure 1(c). 




FIG. 2. H — T phase diagram corresponding to a lattice 
of size L = 12. The vertical straight lino at Tc(L) = 0.84 
corresponds to the L— dependent critical temperature [17], 
which separates the low-temperature ordered phase from the 
high-temperature disordered phase. Open (filled) circles refer 
to the transition between non-defined and concave (convex) 
growth regimes, and squares stand for the Ising-like localizer 
tion-delocalization transition curve. Eight different regions 
are distinguished, as indicated in the figure. Also indicated 
are seven representative points that are discussed in the text. 
The inset shows the phase diagram corresponding to the ther- 
modynamic limit composed of six different regions. 

Further increasing the temperature and for large 
enough fields, the formation of a stable longitudinal up- 
down interface that pushes back the growing interface 
is observed. So, the system adopts the concave growth 
regime (see figure 1(c) corresponding to Region IV in 
figure 2). Increasing the temperature beyond Tc{L), a 
transition from a low-temperature ordered state (Region 
IV) to a high-tcmpcraturc disordered state (Region VI, 
see figure 1(e)), both within the concave growth regime, 
is observed. Analogously, for small enough fields, a tem- 
perature increase drives the system from the ordered con- 
vex growth regime (Region /) to the disordered convex 
growth regime (Region V, see figure 1(d)). As shown in 
figure 2, there is also an intermediate fluctuating state 
(Region B) between Regions V and VI, characterized by 
the competition between the disordered convex growth 
regime and the disordered concave one. 

Finally, a quite unstable and small region (Region ///) 
that exhibits the interplay among the growth regimes of 
the contiguous regions, can also be identified. Since the 
width of Region /// is of the order of the rounding ob- 
served in Tc{L), large fluctuations between ordered and 



disordered states arc observed, as well as from growth 
regimes of non-defined curvature to convex ones. How- 
ever, figure 1(f) shows a snapshot configuration that is 
the fingerprint of Region ///, that may prevail in the 
thermodynamic limit, namely a well defined spin up- 
down interface with an almost fiat growing interface. 

Let us now locate the transition curves between the 
already discussed regions. For this purpose, the mean 
transverse magnetization m{k, L, T, H) is defined as 



1 ^ 

m {k,L,T,H) = j^ ^iih, 



(2) 



and the susceptibility (x) is defined in terms of the mag- 
netization fluctuations. Then, using a standard pro- 
cedure [7], the localization-delocalization "transition" 
curve corresponding to the up-down interface running 
along the walls can be obtained considering that on the 
H — T plane, a point with coordinates {II^,T^) on this 
curve maximizes x{H,T). The obtained curve for the 
confined system is shown in figure 2. This localization- 
delocalization "transition" is also known as a pre-wctting 
"transition" and actually becomes a true wetting transi- 
tion in the thermodynamic limit. 

Since the MEM is a kinetic growth model, another kind 
of nonequilibrium wetting "transition" associated with 
the curvature of the growing interface can also be identi- 
fied. In fact, varying the surface fields one finds a transi- 
tion between a wet state (that corresponds to a concave 
interface) and a nonwet state (associated with a convex 
interface), as already observed in figure 1. Clearly, two 
different contact angles must be defined in order to locate 
this transition, namely 6d for the angle corresponding to 
the dominant spin cluster, and 9md for the one that cor- 
responds to the non-dominant spin cluster. Both contact 
angles can straightforwardly be determined by measur- 
ing the location of the growth interface averaged over a 
sufficiently long growing time. In this way, three differ- 
ent growth regimes can be distinguished: (i) the con- 
cave growth regime that occurs when the system wets 
the walls on both sides (i.e. for 60,6^0 < §), (H) tbc 
convex growth regime that occurs for 9 0,0^0 > f ) and 
(in) the regime of non-defined curvature that occurs oth- 
erwise. The corresponding transition curves obtained for 
the confined system are shown in figure 2. 

We will now extrapolate our results to show that the 
rich variety of phenomena found in a confined geometry 
is still present in the L 00 limit, leading to the phase 
diagram shown in the inset of figure 2. In order to il- 
lustrate the extrapolation procedure, the following seven 
representative points of the finite-size phase diagram are 
discussed in detail: (i) the points labeled Pi, P* , P2, and 
P2, that correspond to the intersections of the if = 0.6 
line with the various transition curves shown in figure 
2, and (ii) the points labeled P3, P3*, and P4, that refer 
to the intersection point between Regions /, //, ///, and 
A, the minimum of the limiting curve between Regions 
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IV-VI and A-B, and the zero-field transition point, re- 
spectively. 
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FIG. 3. Plots of T versus L ^ corresponding to the points 
Pi, Pi , P2, and Pi, all of them with H = 0.6. The fits to the 
data (solid lines) show that Pi P* (i = 1, 2) for L — > 00. 

Figure 3 shows plots of T versus corresponding 

to the points Pi , P* , P2 , and P^ ■ The results from the 
extrapolations are: Ti = 0.67 ± 0.01, = 0.66 ± 0.01, 
and T2 = 1.30 ± 0.02, T* = 1.29 ± 0.01, pointing out 
that, within error bars, Pj ^ P/ (i = 1, 2) in the L — > 00 
limit. Using the same procedure, the extrapolations of 
P3 and P3* (not shown here) give: iJg = 0.30 ± 0.01, 
HI, =0.31±0.02, andr3 = 0.69±0.01,T3* =0.71±0.03. 
So, one has P3 — > P3 for i — > 00 within error bars. Fi- 
nally, the extrapolation of P4 is = Tc = 0.69 ± 0.01. 

Using the above-mentioned extrapolation procedure, 
the phase diagram in the thermodynamic limit can be 
drawn, as shown in the inset of figure 2. By compari- 
son with the finite-size phase diagram of figure 2, notice 
that, besides the fact that the crossover Regions A and B 
vanish in this limit, it is expected that Region /// may 
remain. Although this (very tiny!) region corresponds 
to a physically well characterized growth regime, since 
one expects that the system in this region may grow in 
an ordered phase with a localized up-down domain inter- 
face and a convex growing interface, statistical errors due 
to large fluctuations close to criticality hinder the exact 
location of this region. 

Comparing the equilibrium wetting phase diagram of 
the Ising model [3,6,7] and that of the MEM, it fol- 
lows that the non-equilibrium nature of the latter in- 
troduces new and rich physical features of interest: the 
nonwet (wet) Ising phase splits out into Regions I and II 
(Regions /// and IV), both within the ordered regime 
(T < Tc) but showing an additional transition in the in- 
terface growth mode. Also, the disordered state of the 
Ising system (T > Tc) splits out into Regions V and VI 
exhibiting a transition in the interface growth mode. The 
interplay between wetting of two interfaces, one of them 
defined between differently orientated domains that runs 
along the film, and the remaining one resulting from the 



growing process, leads to a multicritical wetting point 
close to H^'^ = 0.31 ± 0.01, T^^ = 0.70 ± 0.02. At this 
particular multicritical point both interfaces are essen- 
tially perpendicular to each other (figure 1(f)). 

We hope that the presented results will, on the 
one hand, contribute to the understanding of the rich 
and complex physical phenomena exhibited by the irre- 
versible growth of binary mixtures in confined geome- 
tries, and on the other hand, stimulate further experi- 
mental and theoretical work. 
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